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Differential elliptic flow out to p± ~ 5 GeV/c and particle spectra are calculated using the 
MPC elastic parton cascade model for Au+Au at E cm ~ 130 A GeV. The evolution is computed 
from parton transport theory, followed by hadronization either via independent fragmentation 



or by imposing parton-hadron duality. With pQCD elastic cross sections, very large initial 



gluon densities dN g /dr] > 7000 are required to reproduce the data measured by the STAR 
collaboration. In addition, elliptic flow and the p± spectra are shown to be very sensitive to 
particle subdivision. 

° . 

! 1. Introduction 



Elliptic flow, V2(p±) = (cos(20)) p± , the differential second moment of the azimuthal 
momentum distribution, has been the subject of increasing interest Jl|-|5[, especially since 
the discovery [U at RHIC that V2(p± > 1 GeV) — > 0.2. This sizable high-p^ collective 
effect depends strongly on the dynamics in a heavy ion collision and provides important 
^ ■ information about the density and effective energy loss of partons. 

The simplest theoretical framework to study elliptic flow is ideal hydrodynamics^. For 
RHIC energies, ideal hydrodynamics agrees remarkably well with the measured elliptic 
flow data|| up to transverse momenta ~ 1.5 GeV/c. However, it fails to saturate at high 
p± > 2 GeV as does the data reported by STAR at Quark Matter 2001. 

A theoretical problem with ideal hydrodynamics is that it assumes local equilibrium 
throughout the whole evolution. This idealization is marginal for conditions encountered 
in heavy ion collisions 0. A theoretical framework is required that allows for nonequilib- 
rium dynamics. Covariant Boltzmann transport theory provides a convenient framework 
that depends on the local mean free path X(x) = l/an(x). 

Parton cascade simulations]^,^] show on the other hand, that the initial parton density 
based on HIJING|| is too low to produce the observed elliptic flow unless the pQCD 
cross sections are artificially enhanced by a factor ~ 2 — 3. However, gluon saturation 
models predict up to five times higher initial densities, and these may be dense enough 
to generate the observed collective flow even with pQCD elastic cross sections. In this 
study, we explore the dependence of elliptic flow on the initial density and the elastic gg 
cross section. 
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2 



Calculations based on inelastic parton energy lossQ also predict saturation or decreas- 
ing V2 at high p±. These calculations are only valid for high p±, where collective transverse 
flow from lower-px partons can be neglected. The collective component from low p± is, on 
the other hand, automatically incorporated in parton cascades. Though parton cascades 
lack at present covariant inelastic energy loss, elastic energy loss alone may account for 
the observed high-p^ azimuthal flow pattern as long as the number of elastic collisions is 
large enough ||. 

Forerunners of this study ||,|| computed elliptic flow for partonic systems starting from 
initial conditions expected at RHIC. Here, we extend these in three aspects. We compute 
the p^-differential elliptic flow V2(p±), model hadronization, which enables us to compute 
the final observable hadron flow, and consider realistic nuclear geometry. 



2. Covariant parton transport theory 

We consider here, as in P,|T|, |Tt]| , |rT| ] , the simplest form of Lorentz-covariant Boltzmann 
transport theory in which the on-shell phase space density f(x, p), evolves with an elastic 
2 — > 2 rate. We solve the transport equation using the cascade method, which inher- 
ently violates Lorentz covariance. To ensure Lorentz covariance, we apply the parton 
subdivision technique jlO|JT^| . See Q for details. 

The elastic gluon scattering matrix elements in dense parton systems are modeled by 
a Debye-screened form da/dt = cr (l + p? / ' s)p 2 /(t — p> 2 ) 2 , where \x is the screening mass, 
cr = 9na 2 /2p 2 is the total cross section, which we chose to be independent of energy. 



3. Numerical results 



The initial condition was a longitudinally boost invariant Bjorken tube in local thermal 
equilibrium at T = 700 MeV at proper time r = 0.1 fm/c as by fitting the gluon mini- 
jet transverse momentum spectrum predicted by HIJING|§ (without shadowing and jet 
quenching). The pseudo-rapidity distribution was taken as uniform between \r)\ < 5, while 
the transverse density distribution was assumed to be proportional to the binary collision 
distribution for two Woods-Saxon distributions. For collisions at impact parameter b, the 
transverse binary collision profile is dN g (h) / dr]d 2 x.^ = (Tj e tTA{^± + b/2)T^(x^ — b/2), 
where T^(b) = J dz Pa(V z 2 + b 2 ), in terms of the diffuse nuclear density Pa{t)- From 
HIJING, the pQCD jet cross section normalization at b = 0, <jj e t is fixed such that 
dN(0)/dr] = 210. 

The evolution was performed numerically with 40 and 100 mb isotropic cross sections, 
and with 3, 40 and 100 mb gluonic cross sections (/i = T ). We used particle subdivision 
£ = 100 for impact parameters 0, 2, and 4 fm, while £ = 220, 450, 1100, and 5000, for 
b = 6, 8, 10, and 12 fm. 

Two different hadronization schemes were applied. One is based on local parton-hadron 
duality, where as |f§, we assumed that each gluon gets converted to a pion with equal 
probability for the three isospin states. Hence, fh-(p±) ~ fn-(p±) — kfg(p±)- 

The other hadronization prescription was independent fragmentation. We considered 
only the g — > channel with the NLO fragmentation function from [|l^]. For the scale 
factor s = log(Q 2 )/ log(Qo) we take s = because the initial HIJING gluon distribution 
is already quenched due to initial and final state radiation. Since we do not consider soft 
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physics, we limit our study to hadrons with p± > 2 GeV. 

Due to scaling [!]], differential elliptic flow depends on a and dN g /drj only via £ = 
a o dN g (0)/dr]. On the other hand, the p± spectrum depends on a and dN g /drj separately. 

3.1. Elliptic flow results 

Fig. |l] shows elliptic flow as a function of p± and impact parameter. With increasing 
p±, the minimum-bias elliptic flow increases until p± ~ 1.5 — 2 GeV, where it saturates. 
With increasing impact parameter, elliptic flow first monotonically increases, then mono- 
tonically decreases, showing a maximum at b ~ 8 fm. 
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Figure 1. Elliptie<flow as a function of p± for A«HA.u at y/s = 13CL4 GeV. Left figure shows gluon 
elliptic flow for different impact parameters. Middle and right figures show minimum-bias charged hadron 
V2 for hadronization via local parton- hadron duality (middle) or independent fragmentation (right). 



The minimum-bias flow was computed via v™' b '(p±) = J, 71 " Jn max V2(b,p±) b db, where 
we chose b max = 12 fm. This differs from the definition by STAR, which weights V2(b,p±) 
with dN g /dydp± and hence results in a smaller V2(p±) than our estimate. 

The £ = 8400 mb isotropic curve is almost the same as the £ = 21000 mb gluonic curve, 
showing that not the total, rather the transport cross section is the relevant parameter. 
Thus, we expect that the flow in the £ = 21000 mb isotropic case can also reproduced 
with £ ~ 50000 mb and a gluonic cross section. With the pQCD gg cross section of 3mb, 
this corresponds to an 80 times larger initial density than from HIJING. 

3.2. Particle spectra 

Fig. shows the quenching of the p± spectra due to elastic energy loss. In the dilute 
£ = 630 mb case quenching is negligible, however, for £ = 21000 mb it is an order of 
magnitude at p± > 6 GeV for semicentral collisions. Naturally, for a fixed £, quenching 
decreases as the collision becomes more peripheral because the parton density decreases. 
Hadronization via independent fragmentation results in an additional quenching. 

Fig. H] shows that elliptic flow and the p± spectra are sensitive to particle subdivision. 
To obtain covariant results for these conditions, we needed £ > 200. 
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Figure 2. Final 1 gluon and negative hacTron p± spectra (via indepf 'fragmentation) as a fflriSfion of the 
elastic cross section for Au+Au at = 130A GeV with b = 6 fm (left two plots) and b = 12 fm (right 
two plots). For local parton-hadron duality, hadron spectra are proportional to gluon spectra. 




Figure 3. Gluon elliptic flow and p± spectra for Au+Au at y/s — 130A GeV with b = 8 fm and particle 
subdivisions I = 1, 5, 50, 225, and 450. The elastic gg cross section was 100 mb (p. = T ). 
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